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Abstract 

The hadroproduction of charmed mesons is discussed within the framework 
of the modified Quark-Gluon String model (QGSM) taking into account the 
decays of corresponding S'-wave resonances hke D* . It is shown that the esti- 
mate of the cc-pair contribution to the quark sea of colliding hadrons strongly 
depends on the parametrization of the function of the leading light quark frag- 
mentation into charmed mesons. A description of the existing experimental 
data on inclusive spectra and asymmetries of leading and nonleading D - and 
D* - meson production in np -, pp - and T,p - collisions is obtained. The 
predictions for leading/nonleading asymmetry in and collisions 

are also given. 
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1 Introduction 



Hadropro duct ion of charmed particles is now being investigated in many exper- 
iments at various energies of the different initial hadron beams. One of the most 
striking features of charm hadroproduction is the leading particle effect |jl|]. For 
example, in 7r~ (du) interactions with hadrons and nuclei more D~{dc) than D^[cd) 
are produced at large xp- 

Leading charm production can be quantified by studies of the production asym- 
metries between leading and nonleading charm production. The good measurements 
of the differential cross sections up to high value of Feynman xj? — > 1 gives a possi- 
bility of obtaining the leading/nonleading asymmetry defined as 

A(D' D^) - -(D-)--(D'-) a) 
^ - a{D-) + a{D+) 

The experimental data WA82 0, E769 0, E791 § and WA92 § show that the 
asymmetry increases from nearly zero at small xp to 0.5 around xp = 0.6 and does 
not show strong energy dependence in the range pi = 250 — 500GeV/c. This means 
that leading charm asymmetry is primerily located at large xp. Thus the asymmetry 
A reflects the physics at only a small fraction of the total cross section. 

Neutral D's were usually not used in the analysis since they can also be pro- 
duced by the decay of nonleading D*^ - mesons. Recently the WA92 (Beatrice) 
collaboration presented the data on 74(D°,Z)°) asymmetry at 350 GeV/c. 

How can one explain the origin of leading charm asymmetry within the context 
of different theoretical approaches? 

Perturbative QCD at leading order predicts that c and c quarks are produced 
with the same distributions. The asymmetry in this case is equal to zero. 

In PYTHIA, a Monte Carlo program based on the Lund string fragmentation 
model 1^, , it is assumed that heavy quarks are produced in the initial state with a 
relatively small longitudal momentum fraction by the leading twist fusion processes. 
In order to produce a strong leading particle effect at large xp, the string has to 
accelerate the heavy quark as it fragments and to form the final heavy hadron. 

In papers ||, H, [1^, |1T| a QCD mechanism which produced a strong leading/non- 
leading asymmetry at large xp is discussed. The basic assumption of the intrinsic 
charm model |T^, |TT|] is the coalescence of charmed quarks and a projectile valence 
quarks occur in the initial state. For example, vr" can fluctuate into a | udcc > 
Fock state. These states are produced in QCD from amplitudes involving two or 
more gluons, attached to the charmed quarks. Since the charm and valence quarks 
have the same rapidity in these states, it is easy for them to coalesce into charmed 
hadrons and produce leading particle correlation at large xp where this mechanism 
can dominate the production rate. 



The leading effect was also considered in different recombination models |12, 13 



IJ, , where in addition to processes of c-quark hadronization, the recombination 



of the light valence quark with the c-quark is taken into account. 
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In the framework of the Quark-Gluon String model (QGSM) |[T^, |T^, the 
"intrinsic" cc pairs can be taken into account only as a small admixture in the 
quark sea of colliding particles. This assumption is in some disagreement with the 
IC hypothesis ^ |10|, 0] because c and c quarks should rather be considered as 
valence quarks. 

The first attempt to introduce the charmed sea contribution into the QGSM was 
made in |]I9| without comparison with the experimental data on leading/ nonleading 
asymmetry. 

The qualitative analysis of the charmed sea contribution in the framework of the 
QGSM was made in where the term "topological charm" for this phenomenon 
was proposed in order to distinguish it from "intrinsic charm" suggested by Brodsky 
and Vogt 1,0,0]. 

The QGSM prediction for the inclusive spectra of charmed hadrons was consid- 
ered (without intrinsic charm) in |21| - |]26[ . 

In the present paper we analysed the experimental data on leading/nonleading 
D-meson production asymmetry and inclusive cross section in the framefork of the 
modification of the QGSM, developed in [^, ^ taking into account the charmed 
sea contribution. 



2 Model description 

The modification of the QGSM taking into account the contributions from decays 
of corresponding S-wave resonances was developed in p5[. Under this assumption 



the invariant cross section of the produced hadron h integrated over transverse 
momenta p± can be written as 



da da f do- ^. . , . . 

" x*_ 

Here, x '^"^^"^ is the direct production cross section of the hadron h and xr^^ is 
the i?-resonance production cross section. The function ^{xr) describes the decay 
of the resonance R into hadron h. S-wave charmed resonances decay into stable 
charmed particles emitting a vr-meson or a 7-quantum [^, and we describe the 



two-body kinematics of this decay as in [28|. After integration over the transverse 



momenta of both the hadron h and the resonance R the function has the 

form 

In eqs. (|^) and xr is the Feynman variable of the resonance i?. 
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Mrx , Mrx ^ I , 2^<pi>+m2 



+ ~ 171* *' - ~ ZTi* I *' J- — Y "t- "T J^^, _ , 

^ ±j* — p* V 

m is the mass of the produced hadron Mr is the mass of the resonance, E* 
and p* are energy and 3-momentum of hadron h in the resonance rest frame, <p\> 
is the average transverse momentum squared of the hadron h. 

The inclusive spectra of the hadron h in the framework of the QGSM has the 
form 

x-r- = E ('nis)ip^{x) , (4) 

'^^ n=0 

where cr„(x) is the cross section of n-pomeron shower production and f^i^) 
determines the contribution of the diagram with n cut pomerons. The expressions 
for cr„ (s) and the corresponding parameter values for pp and np collisions are given 
in 0-0. 

The function V'^(x) {n > 1) for vrp interaction can be written in the form p3[-p5 



+2(n-l)/Kx+,n)/L(a:_,n)) (5) 
and for baryon-proton interaction 



if^^{x) = a^(/^^(x+,n)/^(x_,n) + n)/^,(x_, n) + 

2(n-l)/L(x+,n)/L(x_,n)) , (6) 



where x± = ^[^x'^ + x]_ ± x]. 

The functions fl^{x, n){i = q, q, qq, g^ea) in © and (H) describe the contributions 
of the valence/sea quarks, antiquarks and diquarks, respectively and were deter- 
mined by convolution of the corresponding distribution functions Ui{x,n) in the 
colliding hadrons with the function of quark/diquark fragmentation into hadron h 
Glix,n) m: 

fi{x,n) = [ Ui{xi,n)Gi{x/xi)dxi . (7) 

The projectile (target) contribution depends on the variable x+ {x-). 
The functions /j^g„(x+,n) and /^g^(x_,n) in (^) and (H) are parametrized in the 
form 
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1 



2 + 6s + Sc 
1 



x± 2 



■dxi 



+ / U2[xi,n)^ 

+ Os u-s[xi,n) dxi 

Jx± Z 

ri ^ G^ix^/x,) + G^,ixJx,) ^ 
+ dc Uc[xi,n) dxi 

Jx± Z 



(8) 



The parameters 6s ~ 0.2 — 0.3 and 6 stand for strange and charmed quark 
supression in the sea. The value of 6c will be found later from the comparison with 
the experimental data. 

A full list of the quark/diquark distribution functions in the vr-meson, p, E~, 



and S -hyperons, used in this work, is given in [25 



We assume the following parametrization of the charmed sea in different hadrons 



(9) 



7: = -ap(0) + (ap(0)-a^(0)). 



P{x,n) = - a;)7?+2{«-i)(i-"P 



(0)) 



(10) 



= K(o) - 2«jv(o)) + K(o) - «^(o)). 



(0)) 



7r = (ap(0) - 2a^(0)) + (^^(O) - a^(0)) + (ap(0) - a^O)). 



u' 



(12) 



7r = («p(0) - 2a^(0)) + (ap(0) - a^(0)) + 2(a^(0) - a^O)). 
The coefficients Cc in (^ - 11^) are determined by the normalization condition 

u^{x, n)dx = 1 



Further we will assume that the fragmentation functions of quarks and diquarks 
do not depend on the spin of the picked up quark (or diquark) This leads to 
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the equality of the fragmentation functions of the corresponding quarks or diquarks 
into different D and D* mesons. In the present paper we use the parametrization of 
the hght quark and diquark fragmentation function given in [21, ^] but with other 
values of free parameters. 

The nonleading fragmentation functions of the quark into D mesons have the 
form 



-a^(0)+A+2(l-aB(0)) 



where A=2a^,(0)pi^., a^(0) = -2.18 [H). 

The leading type fragmentation function is parametrized as: 



(13) 



G 



z)-"'^(°)+\l + arr 



D \ 
1 



(14) 



The aiz^ term was introduced in 
into the /T-meson 



IjJ by analogy with the fragmentation function 
and stands for the density of mesons at the end of the 
fragmenting quark-gluon string. The value of ai is poorly determined from the 



inclusive cross section data (see 



In our approach it is possible to have 



reasonable agreement with the asymmetry data without introducing the charmed 
sea contribution with oi = 2. A bit better agreement can be reached by introducing 
a small fraction of charmed sea( see the next section). Introduction of the term 
1 + 2;^ in ( |14|) is rather arbitrary. We also consider another parametrization 



G^~{z) = (1 - ^)-"^(0)+A(^ _^ ^D^Ay 



(15) 



The fragmentation function of the c(c) - quark into charmed mesons was parametrized 
in the form 



l-a^(O)/ 



\A-aj,(0) 



(16) 



where = 1. 

The value of 6c is different for different parametrizations and was obtained from 
the comparison with the experimental data - - p5| . 

For the hyperon beam processes it is needed to know the fragmentation func- 



tions of the strange quark and the corresponding diquarks. According to |^9[ the 
parametrization was chosen as 



A-«^(0)-afl(0)-a^(0)+2 



(17) 



G^~{z) = (1 - 2)-^+^""'^(°)+°«(°)-2°Jv(0)(- 



l + aiz^ il-z) 



Gg^iz) = (1 - ^)A-a^{0)+2(l-a^(0)^- ^ 



1 (1 - ^) 



(18) 
(19) 
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= (1 - ^)A-a^(0)+2(l-a^(0)) (20) 

The values of the parameters will be given in the next section. 

3 Comparison with the experimental data and 
the predictions of the model 

In this section, we consider the description of the existing experimental data for 
the D and D*-meson hadroproproduction in the framework of the present model 



using fragmentation functions (^41) and (|15|). In what follows, we will present the 
three curves in all figures. They correspond to three versions under consideration: 
1) (i-quark fragmentation function into Z)~-meson ([I^) with the parameters = 
0.0007, af = 2, without taking into account intrinsic charm {5c = 0.). This curve 
is shown by full line and marked by number 1; 2) the same as in 1) but 6c = 0.005, 
dashed line, marked by 2; 3) d-quark fragmentation function into the D~-meson 
according (0) with the parameters = 0.0007, af = 8, 5c = 0.01, dash-dotted 
line, number 3. All sets of parameters were determined from the comparision with 
the experimental data [@ - and |3^ - [0. All theoretical curves in the model 



under consideration are sums of the directly produced Z)-meson cross section and 
the contribution of the decay of the corresponding D* resonanse. 

The experimental data on the xp dependence of {D~ ,D~^) asymmetry measured 
by different groups, WA82 for Pl = 250 GeV/c 0, E769 for Pl = 340 GeV/c [|, 
and the recent data of WA92 (Beatrice collaboration) for Pl = 350 GeV/c are 
presented in Fig.l. The theoretical curves were calculated for Pl = 340 /c. 

The Beatrice collaboration presented the measurements of {D^,D^) asymmetry 

. This data are compared with our calculations in Fig. 2. In Fig. 3 we plotted the 
experimental points for the asymmetry {D~ ,D~^) at Pl = 500 GeV/c together 
with the QGSM calculations in the same kinematic region. 

The experimental data on the xp- dependence of the inclusive distributions of 



all D-mesons in 7r~p and pp interactions at Pl = 200 GeV ^2| and 250 GeV/ c 
are presented in Figs. 4 and 5 respectively. The theoretical curves are calculated for 
the sums of the spectra of all D-mesons at Pl = 250GeV/ c. 

The model calculations for the spectra of the sum of all D mesons in the reaction 
TT^p DX at Pl = 360 GeV/c are compared with the experimental data at 
P^ = 350 GeV/c § and Pl = 360 GeV/c [§D|, §T) in Fig.6. 

As far as we consider here the QGSM spectra of resonances taking into account 
their subsequent decays, it is important to compare our calculations with the avail- 
able data on D*- meson production. The data on the reactions 7[~p D*^ / D*~ X 
and TT^p D*^ / D*^X at 360 GeV/c |^ are compared with our predictions in Fig. 7 



and 8 respectively. 

In they presented the asymmetry of the leading {D*^,D*^) and nonleading 



{D*^, D ) vector resonances with open charm in vr p- collisions, integrated in the 
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xp > region. This value is A{L,NL) = 0.09 ± 0.06. The calculations of the 
model give the following values for all parametrizations used in this work: 1 - 
A{L,NL) = 0.103, 2 - A{L,NL) = 0.097, 3 - A{L,NL) = 0.088. As we can see 
the model calculation agrees with the experimental measurements. The data on the 
xp dependence of differential cross section of the reactions n~p D*~^/D*~X and 
7i~p — i> D*^/D*^X at 360 Gey/c[^] are compared with our calculations in Fig.7 
and 8. These results confirm our assumption about equality of the spectra of the 
vector {D*) and directly produced pseudoscalar {D) mesons. 

The predictions of the model for the inclusive spectra of D~ and mesons in 
Tj^P collision at 330 GeV/c (WA89) are given in Fig.9 and 10. 

The predictions for the xp- dependence of the asymmetry {D~ ,D^) for H^p and 
S'p- collisions are presented in Figs. 10 and 11, respectively. 

As we can see, the existing data do not allow to have definite conclusion on the 
charmed sea contribution in the QGSM. The curves in Figs. 1-3 show that it is pos- 
sible to have a reasonable description of the data on leading/nonleading asymmetry 
without charmed sea contribution if one chooses the value of the corresponding pa- 
rameter ai = 2 (full line in the figures). The charmed sea contribution is noticeable 
in the region xp ^ 1, and negligible for low xp. This is due to large suppression of 
the charmed sea at ^ (x^ °'i'('^'>^^ 

In pO[ the asymmetry behaviour for another set of parameters (ai = 4, 6c = 
O.Olao) was compared with the calculation for the old one which gives 

larger values of the asymmetry for the entire xp interval. 

The comparison of the experimental data with the calculations using parametriza- 
tions of fragmentation function (|T5|) shows that in this case agreement is even slightly 
better. 

The main conclusion from the results considered is that the existing data do not 
allow an unambiguos quantitative answer conserning the charmed sea contribution 
in the QGSM. 

Acknowledgments. The author is grateful to A.B.Kaidalov, K.G.Boreskov, 
E.Chudakov and O.I.Piskounova for useful discussion. This work was supported by 
Grant INTAS-93-0079ext and Armenian Sciense Foundation (Grant 94-681). 
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Figure captions 

Fig.l Comparison of the QGSM calculations with the experimental data on D~ / 

at Pl = 250 GeV/c § and Pl = 340 GeV/c [|] and Pl = 350 GeV/c [|. The 
theoretical curves were calculated for Pl = 340 GeV/c, full curve (1) without 
including charmed sea; dashed (2) with taking into account the charmed sea; 
dash-dotted curve (3) d- quark fragmentation function according the formulae 

(0)- 

Fig. 2 Comparison of the QGSM calculations with the experimental data on D^/D^ 
asymmetry at Pl = 350 GeV/c @. The curves are the same as in fig.l. 

Fig.3 Comparison of the QGSM calculations with the experimental data on D~/D^ 
asymmetry at Pl = 500 GeV/c The curves are the same as in fig.l. 

Fig.4 Differential cross sections of the D-meson production in TT~p- interaction. The 
experimental data are 200 GeV/ (NA32) |2[, 250 GeV/c (E769) [Q. The curves 
are the same as in fig.l. 

Fig. 5 Differential cross sections of the D-meson production in pp- interaction. The 
experimental data are 200 GeV/ (NA32) |2|, 250 GeV/c (E769) [Q. The curves 
are the same as in fig.l. 

Fig. 6 Differential cross sections of the D-meson production in TT~p- interaction. The 
experimental data are 360 GeV/ (NA27) |0|, 350 GeV/c (WA92) §. The 
curves are the same as in fig.l. 

Fig.7 The Xir-dependence of D*^ /D*^ -mesons in tt~p interaction at 360 GeV/c 
The curves are the same as in fig.l. 

Fig.8 The Xi^-dependence of Z}*°/Z)*°-mesons in ir'p interaction at 360 GeV/c ^5 
The curves are the same as in fig.l. 

Fig.9 Model prediction for the Xir-dependence of the D~-meson cross section in T,~p 
interaction at 330 GeV/c. The curves are the same as in fig.l. 

Fig. 10 Model prediction for the XiT'-dependence of the D^-meson cross section in T,~p 
interaction at 330 GeV/c. The curves are the same as in fig.l. 

Fig.ll Predictions for the Xi?-dependence of D~ / asymmetry for the S~ beam at 
330 GeV/c. The curves are the same as in fig.l. 

Fig.l2 Predictions for Xi?-dependence for the / asymmetry on S~ beam at 
330 GeV/c. The curves are the same as in fig.l. 
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